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ABSTRACT 

The m-z relation for Type la supernovae is one of the key pieces of evidence supporting the 
cosmological ‘concordance model’ with Ao ~ 0.7 and Ho ~ 0.3. However, it is well known 
that the m-z relation depends not only on Ao and Ho (with // ( > as a scale factor) but also on the 
density of matter along the line of sight, which is not necessarily the same as the large-scale 
density. I investigate to what extent the measurement of Ao and Ho depends on this density 
when it is characterized by the parameter 77 (0 < r] < 1 ), which describes the ratio of density 
along the line of sight to the overall density. I also discuss what constraints can be placed on 
77 , both with and without constraints on Ao and Ho in addition to those from the m-z relation 
for Type la supernovae. 

Key words: supernovae: general - cosmological parameters - cosmology: theory - dark 
energy - dark matter. 


1 INTRODUCTION 


In the last 15 years or so, cosmological observations have im¬ 
proved greatly and it also appears that the values are converging 
on their true valuesQ This allows us to answer such questions 
(provided, of course, that ‘standard assumptions’ hold) as whether 
the Universe will expand for ever (yes), how old it is, whether it 
is accelerating now (yes), when it started accelerating, etc. (With 
the assumption of a simple topology, the Universe is finite if 
Ao + Ho > 1, but since observations indicate that this value is very 
close to 1, we cannot yet answer this question.) Among the most 
important of these o bservations are those by the Supernova Cos¬ 


mology Project (e . g. Goobar & Perlmutter 


1 19951.119981Il999l: lAmanullah etafhoid : ISuzuki et al 


Perlmutter et al.l 

120121 ) and 


the High-z Supernova Search tea m (e.g. iGarnavich et al.l |l998|; 


Ries s et ~alll 19981: iRiess et alj|2000h (see also the reviews by Riess 
200(1 Leibundsutl 1200 lL ~ 20081 and I Goobar & Leibundgutl 1200 ll) 


which provide joint c onstraints on Ap and Hp. Combined with 
other observations fe.g. lKomatsu et all201 ll : \pianck Collaboration! 
120131 ), these lead to quite well const rained v alues for the cosmo¬ 
logical parameters (e.g. fig. 5 in lSuzuki et alJl20l2h . Although the 
supernova data alone allow a relatively wide range of significantly 


* E-mail: helbig@astro.multivax.de 

1 The case for the ‘conco rdance model’ with A p ~ 0.7 and flo ~ 
0.3 was already made by Ipstriker & Steinhardtl jl995t) : the values of 
the concordance model thus do not need the supernova data, though of 
course adding more data improves the constraints. While th e correspond¬ 
ing uncertainties have dr amatically decreased (e.g. Kom atsu et all l201ll : 
1 Planck Collaborationll2013h . the values themselves have remained constant 
over the last twenty years. 


different other models, it is interesting that the best-fitting values 
obtained from these observations using the current data are quite 
close to the much better constrained values using combinations 
of several observations without the supernova data, at least under 
the assumptions with which the former were calculated. However, 
since the m-z relation depends not only on Ao and Ho (with Ho as 
a scale factor) but also on the distribution of matter along and near 
the line of sight, the dependence of conclusions drawn from the m- 
z relation for Type la supemovae on this matter distribution should 
be investigated. Alternatively, these observations can perhaps tell 
us something about this distribution. 

The plan of this paper is as follows. Section [2] sketches the 
basic theory used in this paper. In section [3] I briefly review pre¬ 
vious investigations of the influence of a locally inhomogeneous 
universe on the m-z relation. Section [4] describes the calculations 
done and discusses the results. Summary, conclusions and outlook 
are presented in section [5] 


2 BASIC THEORY 

iKavser, Helbig & Schramml l l 19971 . hereafter KHS) developed a 
general and practical method for calculating cosmological dis¬ 
tances in the case of a locally inhomogeneous universe. See KHS 
for details (and for a description of the notation, which is followed 
here); here I repeat only the most important points for the purpose 
of this paper. 

If the Universe is homogeneous, then the fact that light prop¬ 
agates along null geodesics provides sufficient information to cal¬ 
culate distances from redshift. If the Universe is locally inhomoge¬ 
neous, then distances which depend on angular observables related 
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to the propagation of radiation will differ from the homogeneous 
case because more or less convergence will change the angle in¬ 
volved (the angle at the observer in the case of the angular-size dis¬ 
tance, that at the source in the case of the luminosity distance). The 
basic idea is that one considers a Universe which is homogeneous 
and isotropic on large scales, this determining the global dynam¬ 
ics via the Friedmann-Lemaitre equation. Local inhomogeneities 
are modelled as clumps, where the extra matter in the clumps is 
taken from the surrounding matter. Thus, a beam which propagates 
between clumps will have only this thinned-out matter inside the 
beam, while outside the beam the average density (taking both the 
thinned-out background matter and the clumps into account) is ap¬ 
proximately equal to the global de nsity ( precisel y so in the limit of 
an infi nites imal beam), Zeldoyic ijj 1964 l. lDashevskii & Zeldovichl 
< 119651) and lDashevskii & Slvshl jl966h developed a general differ¬ 
ential equation for the distance l between two light rays on the 
boundary of a small light cone (the beam) propagating far away 
from all clumps of matter in a locally inhomogeneous universe: 


l = -AnGripl + ^-1 , (1) 

where G is the gravitational constant, R the scale factor, and 77 
(defined below) and the density p are functions of time (a dot in¬ 
dicates differentiation with respect to time). The first term can be 
interpreted as Ricci focusing due to the matter inside the beam, and 
the second term is due to the expansion of space during the light 
propagation. The key assumption here is that while the the density 
p within the beam can differ from the overall density, the over¬ 
all dynamics of the universe is still described by the Friedmann- 
Lemaitre equation. The assumption that the light propagates far 
from all clumps means that Weyl focusing (shear) can be neglected. 
In the case that the densities inside the beam and outside the beam 
are the same, one of course recovers the homogeneous case. 

Since the angular-size distance is defined as D = l / 6 , where 
9 is the angle at the apex of the beam (at the observer, not at the ob¬ 
served object), D follows the same differential equation as l. Mak¬ 
ing use of this, one can derive a general equation for the angular- 
size distance, valid for all (perturbed, in the sense described above) 
Friedmann-Lemaitre cosmological models and all reasonable (see 
below) values of 77 . KHS described the inhomogeneity via the pa¬ 
rameter 0 < 77 < 1 , where 77 is ratio of the density inside the beam 
to the global density or, alternatively, the fraction of matter which 
is homogeneously distributed, as opposed to being clumped^ This 
leads to a second-order differential equation for the angular-size 
distance (equation (33) in KHS) which can be efficiently integrated 
numerically: 

QD " + {jT- + \q') D '+ + = 0 , ( 2 ) 

where 


Q(z) — n 0 (l + z) 3 — (Uq + Aq — 1)(1 + z) T Ao . (3) 


In the locally inhomogeneous case as well the luminosity distance, 
which is needed in this paper, is larger than the angular-size dis¬ 
tance by a factor of (1 + z) 2 . 

This change, compared to the perfectly homogeneous case, is 
essentially a negative gravitational-lensing effect. In a conventional 


2 This is sometimes denoted by a. I. and some others, use 77 because locally 
inhomogeneous cosmological models are often used in gravitational lensing 
(which per se implies local inhomogeneities) where a is almost always used 
to denote the deflection angle. 


gravitational-lensing scenario, if the density at a given redshift be¬ 
tween two light rays is higher than the overall density (the corre¬ 
sponding overdensity being ‘the lens’), then there will be more con¬ 
vergence than in the case where the two densities are the same. In 
the case of light propagating between clumps, as described above, 
the situation is reversed, and the density between the light rays 
defining the distance-related angle is less than the overall density. 
This means that there is (negative) Ricci focusing (and no Weyl 
focusing), making objects appear fainter than they would be in 
the completely homogeneous case. Of course, this is only a rough 
model, but can be expected to be more realistic than the completely 
homogeneous case and to determine not just the sign of the differ¬ 
ence but also give at least an estimate of its strength. 

Obviously, one cannot have 77 < 0. However, it does not make 
sense to have 77 > 1 either. While it is certainly possible that the 
average density inside the beam could be greater than the global 
density, such cases are either unrealistic or not useful. The limit¬ 
ing case where the density in the beam is greater than the global 
density by a constant factor at every redshift is unrealistic because 
this would imply the existence of overdense regions with an ex¬ 
treme length-to-width ratio which are aligned between us and the 
source, which is incompatible with homogeneity and isotropy on 
large scales and would also put us in a special position. The other 
limiting case where a single compact object increases the density 
in the beam to above the global density is certainly possible, but 
observationally would show up as a gravitational-lens effect and 
should be analysed as such (perhaps by adopting 77 ~ 0 for the 
distance calculation and explicitly calculating the amplification). 
Of course, cases between these two extremes are possible, but it is 
clear that 77 must be between 0 and 1 if it is used as an additional 
parameter in the manner described by KHS; lines of sight which, 
due to fluctuations, are slightly denser than the overall density are 
certainly possible, but are not usefully parametriz ed by 77 in the 
style of KHS. (But see lLima. Busti & Santos|[2014l for a toy model 
with an interesting extension of the 77 concept.) Note that 77 does 
not have to be constant as a function of redshift, and the code de¬ 
scribed in KHS supports an arbitrary dependence of 77 on z. Also, 
it could be different for different lines of sight. It was pointed out 
bv lWeinberd il976l) that 77 must be 1 when averaged over all lines 
of sight (allowing for the moment higher-than-global densities to 
be parametrized by 77 > 1 ), which follows from flux conservation. 
However, in practice lines of sight will probably avoid concentra¬ 
tions of matter, due to selection effects or design: distant objects 
will be more difficult to observed if there is luminous matter along 
the line of sight or if there is absorbing matter along the line of 
sight0 If these selection effects do not exist, and if the samp l e is 
large enough, then the ‘Safety in Numbers’ effect illolz & Lindetl 
1200 a) allows one to effectively assume 77 = 1 , with inhomogene- 


3 Matter along the line of sight can increase the apparent brightness and 
thus make objects visible which otherwise would not be. This phenomenon, 
known as ‘amplification bias' in gravitational lensing. is relevant only if the 
luminosity function is steep enough (since otherwise the magnification of 
the area of sky observed, which reduces the number of objects per observed 
area, will dominate, resulting in fewer objects in a flux-limited sample). 
However, the whole point of the m—z relation for Type la supernovae is 
that they are standard candles, or can be adjusted to behave as standard can¬ 
dles with the help of other observations, which means that the differential 
(integral) luminosity function is essentially a delta (Heaviside) function, so 
the amplification bias plays no role here. Also, since objects much fainter 
than supernovae can be detected in the corresponding observations, no re¬ 
alistic amplification would make an otherwise undetectable object visible. 
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ity merely increasing the dispers ion, roughly linearly with redshift. 
However, Idarkson et ali i 2012 l) point out that most narrow-beam 
lines of sight are significantly underdense, even for beams much 
thicker than those considered in this paper. (On the other hand, 
they also point out that this does not necessarily lead to a reduc¬ 
tion in brightness if one drops the assumption that inhomogeneities 
can be modelled as perturbations on a uniformly expanding b ack¬ 
gro und, a point also emphasiz ed bv lBoleiko & Ferreirall2012l : see 
also iBagheri & Schwarzl20l4 ) 

The situation discussed above corresponds to the situation 
where the beam contains a density r/ times the global density at 
a given redshift and outside the beam the density is equal to the 
global density. In practice, this means that a fraction 77 of the mass 
in the Universe is smoothly distributed and a fraction 1 — 77 is con¬ 
tained in clumps outside the beam. Of course, ‘smoothly’ depends 
on the size of the beam; for example, small objects are part of 
the ‘smooth’ component, not only in the limiting case where the 
smooth component consists of free elementary particles. The im¬ 
portant point is that their angular size is small compared to that of 
the beam. 77 is thus also a function of angle: the larger the angle, 
the more representative is the matter within the beam, so that 77 ap¬ 
proaches 1 for large enough angles. Since the beams of supemovae 
at cosmological distances are extremely thin objects (the thinnest 
objects ever studied by science), evidence for 77 < 1 should be 
most obvious in the m-z relation for Type la supernovae. 

A given value of 77 along a given line of sight does not im¬ 
ply that this value does not change along the line of sight, although 
that is of course a possibility, but rather that the influence on angle- 
dependent distances can be described by an effective value of 77 
which is some appropriate average of a value which varies along 
the line of sight. This means that it is possible for the density along 
the line of sight to be larger than the global density at some points, 
but this is not in contrast with the claim above that 77 > 1 is not 
useful as long as the effective value r] e ff < 1. Another complication 
is that essentially all lines of sight to supernovae will have a den¬ 
sity higher than the globally average cosmological density due to 
the overdensities associated with the Milky Way and with the su¬ 
pernova host galaxy (and corresponding clusters)@ However, since 
the absolute magnitudes of supernovae are not known from first 
principles, but rather calibrated from observations, this effect is, to 
a first approximation, unobservable, since it is essentially a renor¬ 
malization of the absolute magnitude. Even if this extra matter as¬ 
sociated with the galaxies at the ends of the beam would increase 
the density inside the beam to larger than the global density, it is 
not useful to think of this as 77 > 1 , since I want to compare the 
standard assumption (completely homogeneous Universe, at least 
as far as light propagation is concerned) with that of a more real¬ 
istic distribution. The point of comparison, the m-z relation for a 
homogeneous Universe, also contains extra matter at each end of 
the beam, and hence extra convergence. As far as I know, no-one 
has ever taken this into account and it is not necessary if one is in¬ 
terested only in the differences. (This would have to be taken into 
account, though, if the absolute magnitude of objects at cosmolog¬ 
ical distances were known independently of observation.) 

Although the term ‘dark matter" suggests something opaque, 
the defining characteristic is lack of interaction with electromag¬ 
netic radiation. Thus, not only does dark matter not glow, it is 
also transparent. It is thus irrelevant whether dark-matter objects 
within the beam significantly cover a source as seen by an ob¬ 


server. (Of course, when comparing observed to calculated bright¬ 
ness, one must correct for extinction due to ‘conventional’ matter - 
it can also be dark in the sense that it does not radiate, but it is 
not transparent.) Here, I am using the term ‘dark matter" to refer 
to the ‘missing matter’, i.e. that responsible for the difference be¬ 
tween the density due to baryonic matter (other non-baryonic but 
known particles (neutrinos) do not increase this significantly) and 
the global density of the Universe, as measured on large scales. Of 
course, non-radiating baryonic matter does exist, but we know from 
constraints from big-bang nucleosynthesis that this cannot be a sig¬ 
nificant fraction of the missing matter. This reflects current usage, 
e.g. the ‘DM’ in ‘ACDM’, and is more convenient than ‘not yet 
identified non-baryonic matter’. 

Since we know that the Universe is not exactly homogeneous 
and isotropic, 77 7 ^ 1 is the most obvious departure from the sim¬ 
plest cosmological model (the Einstein-de Sitter model with Ao = 
0 , Ho = 1 , and 77 = 1 , although the last item is often not stated 
explicitly), but there is not much literature on this topic. (There are, 
though, several recent papers investigating whether ‘dark energy’ 
could be something other than the traditional cosmological con¬ 
stant, e.g. whether the equation of state w differs from — 1 , whether 
it changes with time etc, even though there are no observations 
which indicate this. Of course, that does not mean that one should 
not look.) 

If 77 is allowed to vary from one line of sight to another, 
one could regard this as an additional contribution to the un¬ 
certainty in the distance modulus, much the same as the uncer¬ 
tainty in the absolute magnitude. Theoretically, fitting the obser¬ 
vations for a constant value of 77 would result in a worse fit for 
such cases if this additional uncertainty is ignored or in a larger 
allowed region of parameter space if it is included in the error 
budget. With some assumptions, one could try to take this ad- 
ditional dispersion into account and/ o r correct for it; see, e.g ., 
lAmmullah^Mortsell&Goobarlj 2003^_ GunnarssonrtaL J2006), 
Ijonsson et al.1 d2006l) . Ijonsson, Mortsell & Sollermanl <120091) . In 
practice, with a large number of objects and only a few variables, 
the difference in goodness of fit is well within the expected range 
of values for the case in which 77 is the same along all lines of sight. 
Alternatively, with current data it is also a relatively small contribu¬ 
tion to the error budget. Thus, if observations suggest 0 < 77 < 1, 
it would be unclear if this is evidence for the corresponding value 
of the global value of 77 or whether this is a compromise between 
lines of sight with lower and higher values. However, if observa¬ 
tions indicate 77 = 0 or 77 = 1 , then this would be evidence for the 
corresponding global value, because these are the extreme values 
of 77 and cannot result from averaging. 

Of course, more complicated models are possible. In this pa¬ 
per, I consider only models in which 77 is a constant function of 
redshift and the same along all lines of sigh {3 also, in all cases but 
one, it is independent of the other cosmological parameters. The 
variation between these models, however, is certainly larger than 
the realistic range of the possible influence of 77 on the m-z rela¬ 
tion for Type la supernovae. 


3 BRIEF HISTORY 

The effects of a locally inhomogeneous universe on quantities im¬ 
portant for observational cosmology were first investigated in a se- 


4 I thank Philip Bull for first pointing this out to me. 


5 See iGunnarsson et al.1 (20061) for a discussion of a ^-dependent 77 in the 
context of the m-z relation for Type la supemovae. 
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ries of papers by Zeldovich (1964 

, Dashevskii & Zeldovich (19651) 

and Dashevskii & Slysh (1966). 

Dyer & Roeder (1972) discussed 


the special case of Ao = 0 but with Go as a free parameter for 
p = 0 (wher e ther e is an analytic solution) and for general p val¬ 
ues foyer & Roedej[l973l) . As a result, the distance for p = 0 
is sometimes referred to as the Dyer-Roeder distance. KHS pre¬ 
sented a second-order differential equation and numerical imple¬ 
mentation valid for the general case (—oo < Ao < oo, 0 < 


1969^ Kantowski, Vaughan & Branch 199^_ Kantowski 

1998; 

Kantowski, Kao & Thomas 2000; Kantowski & Thomas 

2001 ; 

Kantowski 2003}) have stressed the importance of p for the inter- 


pretation of the m-z relation for Type la supernovae and have 
provided numerical implementations using elliptic integral^ for 
the special values of p of 0, |, and 1. IPerlmutter et alJ 11999 ) 
considered the effect of p ^ 1 on their results (see their 

fig. 8) and concluded that, at least in the ‘interesting’ region of 
the Aq-Gq parame ter space, it had a negligible effect (see also 
Ijonsson et al.ll2006l ). The reason for the current paper is that, with 
the larger number of supernovae now available, this is no longer 
the case. Further investigation has often been motivated by the m- 
z relation for Type la supernovae (e .g. I Goliath & Mdrtselil200d : 
lM6rtsell.~Go obar & Bergst rdml200lh . It has also been investigated, 
via comparison with explicit ray-tracing through mass distribu¬ 
tions derived from simulations or observations, whether p is a use¬ 
ful parametrization f or local inhomogeneity (e.g. iBergstrom et al.1 
l200d : lMbrtsellll2002 ) (and the conclusion is that it is a useful ap¬ 
proximation, at least for cosmological models which are otherwise 
realistic). 

There seem to be three schools with respect to the attitude 
taken to the possible influence of inhomogeneities on cosmological 
parameters derived from the m-z relation for Type la supernovae. 
One school ignores it completely, assuming a completely homoge¬ 
neous Universe as far as the calculation of the luminosity distance 
is concerned (e.g. lRiess et al.lfl998l) . or provi des some limited jus - 
tification for not considering it further (e.g. IBetoule et al.1120141 ). 
Another school emphasizes that the problem is not completely un¬ 
derstood, the amount of uncertainty's unknown, and even the sign 
of some effects is unclear (e.g. lClarkson et al.l2012l) . A third school 
uses some appro ximation to at least get an id ea of the size of possi- 
ble effects (e.g. lMortsell et all200lh . (While lPerlmutter et alJ 19991 
did consider the possible influence of p, hence belonging to the 
third school, at least at that time, with their data then it was not a 
significant source of uncertainty in their main result. One purpose 
of this paper is to show that this is no longer the case.) 


4 CALCULATIONS, RESULTS AND DISCUSSION 


I have used the publicly available ‘Union2.1’ sample of supernova 
data dSuzuki e t alj201 2h and calculated \ 2 and the associated prob¬ 
ability following Amanullah et alJ d2010h on regularly-spaced grids 
of various extents and resolutions in the Ao-Go-r? parameter space. 
This assumes , of course, that p is a free parameter on the same foot¬ 
ing as Ao and Oo- My goal is not to obtain the ‘best’ cosmological 
parameters, not even the ‘best" ones from the supernova data alone. 
Rather, it is to investigate the influence of p ^ 1 on the interpreta¬ 
tion of the m-z relation for Type la supemovae. I have thus inten¬ 
tionally made the supernova data as precise as possible, by using 
only the statistical uncertainties (i.e. column 4 in the publicly avail¬ 
able data file) and fixing Hq at 70 kms _1 Mpc _1 , which implies 
M = —19.3182761161. Thus, all increase in the allowed region 


of parameter space (at a given confidence level) is due only to the 
influence of 7)|j 

I have calculated y 2 and the corresponding probability on two 
three-dimensional grids: a larger, lower-resolution grid 


-5< Ao < 5 AA 0 =0.02 
0 < U 0 < 10 AGo = 0.02 
0 < p < 1 Ap =0.01 


(500 values) 
(500 values) 
(100 values) 


and a smaller, higher-resolution grid 


0 < A 0 < 1.5 AA 0 =0.003125 
0 < Go < 1-0 AGo =0.003125 
0 < p <1.0 A p = 0.01 


(480 values) 
(320 values) 
(100 values) 


(This paper contains no plots based on the larger, lower-resolution 
grid; the corresponding calculations were done to make sure that 
there is no appreciable probability outside of the range of the 
smaller, higher-resolution grid.) Since three-dimensional contours 
cannot be fully represented in two dimensions, I present various 
two- and one-dimensional visualizations in order to illustrate the 
influence of p. 

All contours in two (three) dimensions have been calculated 
as the smallest-area closed curve (smallest-volume closed surface) 
which encloses the corresponding fraction of the probability. I have 
used the standard values 0.683, 0.954 and 0.997; these correspond 
to 1—er, 2—cr and 3-cr in the Gaussian case. However, I have made 
no assumption about Gaussianity, since I have calculated the con¬ 
tours explicitly, rather than plotting them at the corresponding frac¬ 
tion of the peak likelihood under the Gaussian assumption. For all 
plots, the area outside of the plot has been assigned a probability 
of zero. Otherwise, no priors other than those explicitly stated have 
been used. In particular, no prior information on the values of the 
cosmological parameters from other tests have been used; what I 
show depends on the supernova data only. 

Figs. □ 0 and [3] show projections of the three-dimensional 
contours along one axis on to the plane spanned by the other two 
axes for the smaller, higher-resolution grid. It can be seen that the 
combination of Ao and Go is well constrained, as are both individ¬ 
ually, while r] is hardly constrained at all. Note also that Ao and Go 
are less constrained for lower values of p. (The contours at 0.954 
and 0.997 cannot be distinguished in these plots.) The relatively 
sharp bend in the lower right contours in Figs.Q]and[2]is due to the 
fact that I have assigned a probability of 0 to modelswhich have no 
big bang (see the discussion of figs. 1 and 2 in iHelbisj d2012h and 
references therein for an explanation). 

Another way of visualizing these three-dimensional contours 
is to make cuts through them for a fixed value of one of the param¬ 
eters. Figs.[4][5] and[5]show cuts for p = 0.005, 0.455, 0.955. The 
contours become smaller and move to lower values of Ao and Go 
as p becomes larger. (Again, the contours at 0.954 and 0.997 can¬ 
not be distinguished in these plots.) Fig. [6] is quite similar to stan¬ 
dard presentations of the supernova constraints (e.g. ISuzuki et al.1 


6 Since the goal is not to obtain the best constraints on Ao and Go, but 
rather to investigate the influence of p on the constraints. I have retained 
the Union 2.1 sample with which I began this inv estiga tion, rather than up¬ 
dating it to use, e.g., that used by IBetoule et all {2014). Those with better 
access to such data will always h ave a b etter sample than that which is pub¬ 
licly available. Since even IBetoule et al] 1 2014) do not consider p at all. it 
is perhaps important at the moment for a theorist to take a step back for a 
more general view in order to contrast with continual updates using some¬ 
what better samples. It is important, though, that the Union 2.1 sample is 
significantly larger than those used in the early works discussed in section|T| 
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x 


Figure 1. Projection of three-dimensional probability distribution along the 
r/-axis. 



X 


Figure 2. Projection of three-dimensional probability distribution along the 

f2o _ax i s - 

l2012h . but keep in mind that these contours are a cut through 
the three-dimensional contours for a fixed value of 77 , not two- 
dimensional contours. If 77 is substantially less than 1, then not only 
is the allowed region much larger, but the ‘concordance model’ 
with Ao « 0.7 and f2o ~ 0.3 is ruled out. Qualitatively, this be¬ 
haviour is easy to understand: there is some degeneracy between 
77 and Ao + flo since both increase the amount of focusing in the 
beam, the former because there is more matter in the beam and the 
latter because of the increase in the global curvature, which is es¬ 
sentially Ao + flo- When there is essentially no matter in the beam, 



a 

Figure 3. Projection of three-dimensional probability distribution along the 
Ao-axis. 



X 


Figure 4. Cut through the three-dimensional probability distribution per¬ 
pendicular to the 77-axis for r] = 0.005. 

then the value of flo is less important and hence not as well con¬ 
strained. This means that Ao + fio can be realized via a larger range 
of each parameter, making the allowed region larger. The middle 
value of 77 is that of the global maximum probability. (Since Ao and 
S7o are better constrained, the corresponding plots for fixed values 
of these parameters, not shown here, are less interesting.) 

The 'standard procedure’ for reducing the number of param¬ 
eters shown in a plot is to marginalize over the less interest¬ 
ing or ‘nuisance’ parameters. This is shown in Figs. 13 ID and [3 
Here, and in similar figures below, the grey-scale corresponds to 
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x 


Figure 5. Cut through the three-dimensional probability distribution per¬ 
pendicular to the //-axis for r/ = 0.455. 



X 


Figure 6. Cut through the three-dimensional probability distribution per¬ 
pendicular to the p-axis for r/ = 0.955. 

the probability^ These are qualitatively similar to the projections. 
Fig. [7] also contains two straight lines corresponding to a flat uni- 

' It has become fashionable to plot contours and have the regions between 
the contours filled with a certain colour (or perhaps shade of grey). This 
conveys no information in addition to the contours themselves. Of course, 
the probability between two contours, or within the smallest contour, is not 
everywhere the same, as is obvious from Fig. [7] I have chosen to display 
this potentially important information in addition to the contour curves. 



X 


Figure 7. Two-dimensional probability distribution obtained by marginal¬ 
izing over r\. 
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X 

Figure 8. Two-dimensional probability distribution obtained by marginal¬ 
izing over Q 0 . 


verse with Ao + flo = 1 (negative slope) and zero acceleration 
(go — — Ao = 0) (positive slope). Note that a flat universe 

is compatible with the data but not required by them; in fact, the 
degeneracy in the constraints is almost perpendicular to the flat- 
universe line. In this particular plot, the degeneracy corresponds 
roughly to go « —0.6; in many of the other plots, the degeneracy 
in the Ao-Oo plane is closer to a constant value of fio — Ao than 
to a constant value of go = — Ao. (go was important histori¬ 

cally since the departure from the linearity of the m-z relation at 
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Figure 9. Two-dimensional probability distribution obtained by marginal¬ 
izing over Aq . 



X 


Figure 10. Two-dimensional probability distribution obtained by maximiz¬ 
ing T], 


low redshift is proportional to qo \ nowadays quoting a value for qo 
derived from the m-z relation for higher-redshift objects is neither 
necessary nor sufficient nor, in general, meaningful.) 

Another approach is to maximize the ‘nuisance’ parameter, 
i.e. for a given point in the plane of the plot, find the value of the 
third parameter which maximizes the probability. This is shown in 
Fig-llQl (For these data, such plots are very similar to those where 
the third parameter has been marginalized over, so only this one 
example is shown.) 



X 

Figure 11. Two-dimensional probability distribution for q = 0. 



X 

Figure 12. Two-dimensional probability distribution for r/ = 0.455. 

Most discussion of the m-z relation for Type la supernovae 
has concentrated not on contours of more than two dimensions, 
nor on some reduction (projection, cut, marginalization, maximiza¬ 
tion) of these higher-dimensional contours to two dimensions, but 
rather on two-dimensional contours, i.e. with a (5-function prior 
on the nuisance parameters. Almost always, of course, the (often 
implicitly assumed) prior is q = 1. For comparison, in Figs. IIII 
ED and uni show constraints in the Ao-f2o plane for fixed val¬ 
ues of q, namely 0, 0.455 (the value at the maximum of the three- 
dimensional probability distri butio n) and 1. The last should be 
compared with e.g. fig. 11 in [Kowalski et alj J2008I ). but keep in 
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x 

Figure 13. Two-dimensional probability distribution for 77 = 1. 

mind that, as mentioned above, I have fixed Ho and useonly the sta¬ 
tistical uncertainties. (See also figs, la and 5a in lAmanullah et all 
fc003lU Thus, Fig. ED has slightly smaller contours than similar 
plots elsewhere in the literature. Again, this is intentional so that 
any deviations from this fiducial plot (larger and/or shifted con¬ 
tours) are due solely to the influence of 77 . 

Of course, little significance should be placed on variations 
in the probability within the innermost contour, since the proba¬ 
bility that the point representing the true values of Ao and fio is 
only about twice as likely to lie inside this contour than outside it. 
Nevertheless, it is remarkable that the maximum of the probabil¬ 
ity in Fig. |T3]is at Ao = 0.7210938 and Ho = 0.2773438, i.e. at 
the values of the concordance model (within the small uncertain¬ 
ties; these are much smaller than even the 68.3 per cent contour 
in Fig.ll3tFlNote that when fewer supernova data were available, 
the best-fitting value was a t much higher values of Ao and Ho; see 
e.g. fig. 1 in lllelbii] Il999l ). (As mentioned above, the best-fitting 
value is often not visible in modern versions of such plots, though 
of course it can be easily found in the data used to make the plots.) 
If the best-fitting value remains the same when significantly more 
supernova data are available, then very probably the true value will 
have been converged upon, even though the range of values al¬ 
lowed, even at the 68.3 per cent level, would include values well 
outside what is acceptable when other cosmological constraints are 
considered (i.e. joint constraints from several cosmological tests). 
Normally, when more data are available one expects the new best¬ 
fitting value to be consistent with, but different from, the old best¬ 
fitting value, as has been the case with the supernova data up until 
now. However, looking towards the future, I don’t expect the best¬ 
fitting values for Ao and Ho to change significantly, but do expect 
the constraints from the supernova data to improve, which appears 
somewhat puzzling. A possible explanation for this is that the sta- 

8 For completeness, I quote the exact position of the maximum as calcu¬ 
lated on the grid; of course, this does not imply that the maximum is known 
to greater precision than the resolution of the grid. 



X 


Figure 14. Two-dimensional probability distribution for k = 0. 


tistical errors in the supernova data have been overestimated. Note, 
however, that the best-fitting values for the supernova data corre¬ 
spond to the concordance model only if one assumes 77 ~ 1. For 
77 = 0.455, the concordance model lies very near the 95.4 per cent 
contour, and for 77 = 0 it is even outside the 99.7 per cent contour. 
The plots above illustrate that it is not possible to appreciably con¬ 
strain 77 from the supernova data alone. However, the fact that the 
supernova data suggest the concordance model only for high values 
of 77 could be seen as evidence that 77 « 1 . 

A similar result is shown in Fig. ED where a flat universe (Ao + 
Ho = 1 ) has been assumed. As in the other plots, Ao is reasonably 
well constrained, while 77 is quite weakly constrained. (In this case, 
since Ho = 1 — Ao, Ho is just as well constrained; in general, Ho 
is less well constrained than Ao.) However, note that the best-fitting 
value is for 77 = 1 and Ao sa 0.72; in other words, again the best fit 
is for the concordance model with 77 = 1. (This plot also shows the 
importance of plotting the probability and not just a few contours.) 

To illustrate the change in the effect of 77 now that more su¬ 
pernova data are available, Fig.ll5lshows the constraints where 77 is 
a function of Ho. namely 77 = 0 for Ho < 0.25 and 1 — 0.25 f or 
Ho > 0.25. This should be compared with fig. 8 in lPerlmutter et alj 
Jl999l) . In that figure, the red contours were calculated in the same 
way as those in Fig. ED In the same figure, the green contours were 
calculated in the same way as in Fig. QT| The comparison illus¬ 
trates vividly the fact that the effect of 77 can no longer be neglected. 
While [Perlmutter et alj 11999h concluded that, at least in the inter¬ 
esting part of parameter space, the constraints on Ao and Ho from 
the supernova data did not depend heavily on the assumed value of 
77 , this is definitely no longer the case. 

While the supernova data cannot usefully constrain 77 , as has 
been shown above, the fact that they result in the concordance 
model if one assumes 77 ~ 1 suggests that 77 « 1. Since there 
are many cosmological tests completely independent of the super¬ 
nova data, and also independent of the value of 77 , which suggest 
the concordance model (this is of course why it is called the con¬ 
cordance model), one can assume the concordance values for Ao 
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Figure 15. Two-dimensional probability distribution with 77 = /(f2o)- 
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Figure 16. One-dimensional probability distribution for the concordance 
model. 

and and calculate the probability of 77 from the supernova data 
with these additional constraints; this is shown in Fig. ri 6 l The best- 
fitting value is 77 = 0.7485 while the formal statistical limits are: 

0.60 < 7] < 0.90 (68.3%) 

0.46 < 7] < 1.00 (95.4%) 

0.28 < 7] < 1.00 (99.7%) 

. (This can be contrasted with Fig. DU which shows the value 
of 77 preferred by the supernova data alone; Ao and Do have been 
marginalized over.) While 77 = 1 is not ruled out at high confidence, 



Figure 17. One-dimensional probability distribution after marginalizing 
over Ao and Qo ■ 

lower values of 77 are ruled out at a high level of statistical signif- 
icance[f| This suggests that 77 is relatively large, even though the 
beams of supernovae at cosmological distances are extremely thin 
and this cosmological test should suggest a value of 77 lower than 
that of any other cosmological test known today. One would not 
expect to obtain 77 = 1 since some matter is associated with galax¬ 
ies which are outside the beam; such mass contributes about 0.1 to 
Do. Fig.[l 6 ]thus suggests that dark matter is distributed much more 
smoothly than galaxies. While the beam of a supernova at cosmo¬ 
logical distance is almost a fair sample of the Universe, it is an even 
fairer sample of dark matter. Dark matter is thus not significantly 
clumped at the scale of a supernova beam. 

5 SUMMARY, CONCLUSIONS AND OUTLOOK 

The following conclusions were more or less expected. 

(i) Constraints on Ao and Do are weaker if 77 is not constrained. 

(ii) The concordance model is reasonably probable. 

(iii) There is a degeneracy between 77 and the amount of spatial 
curvature (Ao + Do). 

(iv) Ao is constrained best, then Do, then 77 . 

The following conclusion was neither expected nor surprising. 

9 This can be contrasted with the location of the maximum of the three- 
dimensional probability distribution, where the best-fitting values are Ao = 
0.8609375, D 0 = 0.5015625, and 77 = 0.455. While this lies outside 
the allowed region of parameter space as determined from cosmological 
tests other than the m-z relation for Type la supemovae, the allowed re¬ 
gion is quite large and the concordance model with 77 = 1 is within the 
68.3 per cent contour. Even though the constraints on Ao and Do are of 
course weaker if 77 is allowed to vary, a significant portion of the three- 
dimensional parameter space can be ruled out. and portions of the Ao—Do 
plane are also ruled out, though no additional region is ruled out which is 
not already ruled out by other cosmological tests. 
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(i) Even when 77 is allowed to be a free parameter, the m-z rela¬ 
tion for Type la supernovae is not compatible with go = — Ao > 
0, and thus implies that the Universe is currently accelerating^ 
(Even though the m-z relation for Type la supernovae is one of the 
key pieces of evidence supporting the cosmological ‘concordance 
model’ with Ao ~ 0.7 and Ho ~ 0.3, it is not an essential piece 
in the sense that combinations of other tests still result in the same 
concordance model. Nevertheless, it is still an important piece of 
evidence in favour of the concordance model since it is the only 
single test which, without additional assumptions, implies go < 0, 
i.e. a Universe which is currently accelerating.) 

The following conclusions are somewhat surprising. 

(i) The overall (in the three-dimensional parameter space) best¬ 
fitting values for Ao and Ho are ruled out by other cosmological 
tests. Probably, this best-fitting point is the result of overfitting: its 
probability is not significantly higher than elsewhere and the al¬ 
lowed region is quite large. 

(ii) If one assumes k = 0 then the best fit is very close to the 
concordance model and has 77 = 1. 

(iii) If one assumes 77 = 1, then the best fit is very close to the 
concordance model. 

(iv) If one assumes the concordance model, then one can probably 
rule out low values of 77 , even though the relevant scale is extremely 
small, which implies that dark matter is much less clustered than 
galaxies are. 

(v) We cannot rule out 77 = 1, and there is some tentative evidence 
for it. 

To summarize, allowing 77 , which is otherwise only weakly 
constrained, as a free parameter significantly alters both the best 
fit in the Ao-Ho plane and the allowed region of this plane. The 
concordance model is, however, still allowed. There are hints that 
77 ss 1 , though these are not statistically significant when exam¬ 
ined in the three- or two-dimensional parameter space. On the other 
hand, if one assumes the concordance values for Ao and Ho, low 
values of 77 can probably be ruled out. which is not obvious consid¬ 
ering the very small scales involved; this implies that dark matter is 
very homogeneously distributed. 

One might h ave thought that the i ncrease in the number of 
data points since IPerlmutter et alj i ll 9991) would allow some sort 
of useful constraint to be placed on 77 from the supernova data 
without further assumptions. This is not the case. Even worse, if 
77 is allowed to vary, then the conclusions about the cosmologi¬ 
cal model derived from the m-z relation for Type la supemovae 
are not as robust. However, as discussed in section [T| current con¬ 
straints from combinations of cosmological tests without using the 
supernova data determine the ‘concordance model" with Ao S 3 0.7 
and Ho ~ 0.3 to rather high precision. It is thus perhaps more 
interesting to assume the concordance model and use the super¬ 
nova data to constrain 77 , especially since 77 is otherwise difficult 
to measure. Indeed, as shown in Fig. M current data already pro¬ 
vide interesting constraints. It is also extremely interesting that the 
supernova data have the best-fitting values for Ao and Ho corre¬ 
sponding to those of the concordance model if and only if 77 ~ 1 
is assumed. (Note that while the best-fitting value of 77 assuming 
the concordance model is « 0.75, the best-fitting values of Ao und 


10 Mortsell & Clarkson (2009) have shown that this conclusion also holds 
for a much wider class of models than the Friedmann-Lemaitre models 
considered here. 


Ho assuming 77 R 3 0.75 are different from those of the concor¬ 
dance model.) If this is not a statistical fluke, it could indicate that 
77 « 1 , which is somewhat surprising since the value of 77 as ‘felt’ 
by the supernova might be expected to be somewhat less, because 
the corresponding beams are extremely thin. The fact that even the 
supernova data ‘want’ 77 ~ 1 could indicate that dark matter is dis¬ 
tributed extremely homogeneously. See iHolzl J 19981 ) for a different 
expression of the same idea. Alternatively, this could be evidence 
that the ‘Safety in Numbers’ scenario mentioned in section[2]is in 
fact a valid approximation. 

In contrast to the first useful determinatio ns of Ao and Ho from 
the m-z relation for Type la supernovae e.g. lGarnavich et alj 19981 ; 
iRiess et aTIl 1 9981 : [Perl mutter et al.ll 1999L where the effect of 77 ^ 1 
had a negligible effect on the constraints derived, at least in the 
‘interesting’ region of the Ao-Ho parameter space, with the larger 
number of supernovae now available, this is no longer the case. 
At the same time, current supernova data alone cannot usefully 
constrain 77 (though this might be possible if other cosmological 
data are taken into consideration, as discussed in the previous para¬ 
graph). This should be taken into account in attempts to determine 
further parameters, such as w, the equation-of-state parameter for 
dark energy. When more supernova data become available, espe¬ 
cially at higher redshift, it might be possible to usefully constrain 
77 and/or discriminate between the effect of 77 and other parameters 
such as w. (The difference in apparent magnitude for different val¬ 
ues of 77 increases with increasing redshift, while the difference due 
to different values of Ao and Ho is stronger (than that due to varia¬ 
tion in 77 ) at lower redshift and, for some sets of models, decreases 
at higher redshift.) While allowing 77 to be a free parameter, but 
constant as a function of redshift and for different lines of sight, 
is certainly not the last word with respect to the influence of lo¬ 
cally inhomogeneous cosmological models on the m-z relation for 
Type la supernovae, it does demonstrate that care is needed when 
interpreting conclusions derived from assuming 77 = 1. At least, the 
uncertainty in Ao and Ho must be correspondingly increased. While 
it might be possible to decrease this with a more realistic model, it 
is no longer possible to assume 77 = 1 and have confidence in the 
parameters and their uncertainties resulting from an analysis of the 
m-z relation for Type la supemovae. 
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